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Abstract 

This paper presents an experimental, numerical and analytical study of the effect 
of open roof on acoustic propagation along a 3D urban canyon. The experimental 
study is led by means of a street scale model. The numerical results are performed 
with a 2D Finite Difference in Time Domain approach adapted to take into account 
the acoustic radiation losses due to the street open roof. An analytical model, based 
on the modal decomposition of the pressure field in a horizontal plane mixed with 
a 2D image sources model to describe the attenuation along the street, is also 
proposed. Results are given for several frequencies in the low frequency domain 
(1000 — 2500 Hz). The comparison of the three approaches shows a good agreement 
until / = 100 Hz at full scale, the analytical model and the 2D numerical simulation 
adapted to 3D permit to modelize the acoustic propagation along a street. For higher 
frequency, experimental results show that the leakeage, due to the street open roof, 
is not anymore uniformly distributed on all modes of the street. The notion of leaky 
modes must be introduced to modelize the acoustic propagation in a street canyon. 

Key words: urban acoustics, street canyon, Finite Difference in Time Domain 
method, scale model, modal decomposition. 
PACS: 43.50. Vt, 43.28.En, 43.20.Y1, 43.20.E1. 



1 Introduction 



Urban acoustic researches are divided in three thematics : the sources charac- 
terization and identification, the acoustic propagation and the noise perception 
in an urban context. This work enters in the second thematic : its purpose is 
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to describe the propagation of sound emitted by known sources in a street by 
taking into account its physical morphology (open roof, height and width of 
the street). 

To study acoustic propagation in an urban context, several approaches are 
available : energetical methods [1,2,3,4] based on the estimation of a quadratic 
quantity (energy density or acoustic intensity), numerical methods [5,6] to es- 
timate acoustic pressure or velocity and modal approach [7] to calculate the 
pressure or velocity fields. 

Energetical methods are largely used in urban acoustics but concern a limited 
frequency range : the image sources method [1,2,8,9,10], the ray tracing ap- 
proach [11], the radiosity method [12,13,14] and finally statistical methods of 
particle transport [3,4,15] are generally used for middle and high frequencies. 
All these approaches propose to modelize the effect of the street open roof by 
a complete absorption of the wave. 

The numerical methods, as Finite Element Method (FEM) or Boundary El- 
ement Method (BEM), are restricted for urban acoustics to low frequencies 
for 2D case or very low frequencies for 3D case [6,5] because of the large time 
cost dur to the discretization. 

The modal approach, where the geometrical characteristics of the street are 
explicitely taken into account in the model, is generally not used due to the 
complexity of the medium and to the difficulty to determine the modes of an 
open space like a street canyon. For example, Bullen and al [7] have studied 
the acoustic propagation in a guide with infinite height or more recently, the 
modal approach was used to calculate the 2D field in a street section, the 
acoustic radiation conditions being described by an equivalent sources distri- 
bution at the interface between the street and the free space [16]. 

This review highlights more particularly that the open aperture of the street 
on the half free space, essential characteristic of the urban environment, is 
taken into account only in 2D (in a section) and rarely in 3D urban acoustic 
problem [17]. Only the energetical approaches modelize this characteristic by 
a complete absorption of the wave but these methods are restricted to middle 
and high frequencies and this asumption of complete absorption should be 
justified at these frequencies and moreover at low frequencuencies. 

The aim of this paper is to study the acoustic propagation along a 3D street 
canyon and to propose models to describe the role of the street open roof 
in the propagation for low frequencies. This work compares experimentally, 
analytically and numerically the study of the open roof effect on acoustic 
propagation. Experimental results, obtained with a street scale model are 
described in section 2. The analytical study, based on the modal decomposition 
of the field in the horizontal plane mixed with an image sources model to 
determine the attenuation due to the street open roof is presented in section 
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3. The numerical study, made with a 2D-FDTD computation adapted to take 
into account the acoustic radiation losses by the street open roof is proposed in 
section 4. Then, in section 5, experimental, analytical and numerical pressure 
field maps are compared for several frequencies providing validity and limits 
of the models. 



2 Experimental study of the open roof effect on the acoustic prop- 
agation in a street canyon 

2. 1 Experimental set-up 

2.1.1 Scale model of the street 

A 1/25 scale model of 0.27 m height and 0.2 m width with 3 m length is 
carried out (see fig. 1) corresponding in full scale to a 6.75 m heigth, 5 m 
width and 75 m length street. This scale model is fully described in [18] where 
the building facades are simulated by wood cubes with a plane distribution. 
The scale model is put in a semi-anechoic room. 

B&K preamplifier 




On one side, the source is enclosed (see figure 1) in a flat rigid wall while an 
anechoic termination is carried out on the other side (fig. 2b). This anechoic 
termination is made of melamine dihedron designed to obtain a cut-off fre- 
quency around 750 Hz. This allows to consider the scale model of 3 m length 
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as semi infinite for acoustic frequency upper than 1kHz. 




Fig. 2. Views of the street scale model, (a) Rigid wall with a circular aperture at 
the entry of the street scale model, (b) Street scale model with the anechoic end. 

2.1.2 Sound source 

The acoustic source in the street is a loudspeaker enclosed in a rigid box 
opened on a guide with a square cross-section (0.05 x 0.05 m 2 ). Two sorts 
of aperture in the entry rigid wall are used for the experimental studies : a 
0.05 x 0.05 m 2 square cross-section aperture and a circular aperture with 0.01 
m diameter (see the fig. 2a). The square cross-section aperture simulates a 
plane source until 3400 Hz and the circular source simulates a point source. 

2.1.3 Data acquisition system and post-processing 

The acoustic pressure is measured in the speaker box and in the street by 
means of two 1/4 in. pressure microphones (B&K 4938) connected to a pream- 
plifier (B&K 2670) and a conditionning amplifier (B&K Nexus 2693). The 
preamplifier with the microphone is put vertically in the scale model to mini- 
mize the acoustic diffraction at high frequencies and for practical use during 
displacement. The excitation signal is sinusoidal with a variable frequency 
/. The acquisition of the acoustic pressure is performed using a sampling fre- 
quency F s and during a time length T a . The RMS value of the acoustic pressure 
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is determined for each point of the defined scale model mesh. By means of a 
Charlyrobot 3D-robot, a map of the RMS acoustic pressure is then available. 



2.2 Experimental results 

In this work, the acoustic propagation along a street canyon is led for low 
frequency case. The study of the pressure field in a horizontal plane along the 
street is made for the frequencies / = 1000, 1500, 2000, and 2500 Hz. This 
frequencies correspond to / = 40, 60, 80, and 100 Hz in full scale for a street 
of 6.75 m height, 5 m width and 75 m length. The scale model of the street 
is L = 2.8 m length, d = 0.2 m width and h = 0.27 m heigth. The 3D robot 
allows to obtain horizontal maps of the acoustic pressure RMS value along 
the street. The step of spatial sampling is 0.01 m on the x-axis and y-axis. 
The sampling frequency F s is chosen as F s = 20/ and the acquisition time of 
the acoustic signal for each point of the map is defined by T a = N s /F s with 
N s the samples number (typically N s = 2000 which provides 100 samples per 
acoustic period). The RMS value of the acoustic pressure is determined by 
means of a Matlab program using a least mean square method to determine 
the mean value, the amplitude and the phase of the signal. The normalized (to 
the entry) acoustic pressure RMS value (acoustic pressure in the following) 
along the street will be shown. 

The fig. 3a shows the acoustic pressure map along the street for / = 1000 
Hz, with a square source centered on y s = 0.175 m and z s = 0.07 m. The 
map is measured at the height z = 0.07 m. Firstly, the attenuation increases 
along the street. This phenomenon shows that the open roof of the street 
brings acoustic losses for the pressure field inside the street. This attenuation 
is studied more precisely in the following section. Secondly, the shape of the 
acoustic map shows the presence of acoustic modes in the (x, y) plane of the 
street. For / = 1000 Hz, this pressure map can be compared to a map of the 
acoustic propagation in a 2D infinite waveguide with two propagative modes 
without attenuation. 

The fig. 3b shows acoustic pressure along the street at z — 0.07 m for / = 1500 
Hz, with a square source centered on y s = 0.175 m and z s = 0.07 m. The same 
remarks as for / = 1000 Hz can be made, except that the shape of the pressure 
map is not the same due to a different repartition of the source condition on 
the modes. We can note that attenuation depends on frequency showing that 
the acoustic losses due to the open roof depends also on frequency. The fig. 3c 
shows the acoustic pressure map along the street at z = 0.07 m for / = 2000 
Hz, with a point source centered on y s = 0.115 m and z s = 0.07 m. The modal 
repartition of the source condition and the number of the propagative modes 
differ from the two previous cases (for / = 2000 Hz, 3 modes are propagative). 
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Fig. 3. Acoustic pressure in the street (scale model) for (a) / = 1000 Hz at z = 0.07 
m and for a square source centered on y s = 0.175 m and z s = 0.07 m, (b) / = 1500 
Hz at z = 0.07 m and for a square source centered on y s = 0.175 m and z s = 0.07 
m, (c) / = 2000 Hz at z = 0.07 m and for a point source centered on y s = 0.115 cm 
and z s = 0.07 m, (d) / = 2500 Hz at z = 0.07 cm and for a square source centered 
on y s = 0.1 m and z s = 0.07. 



The fig. 3d shows the acoustic pressure attenuation along the street at z = 0.07 
m for / = 2500 Hz, with a point source centered on ?/ s = 0.1 m and z s = 0.07. 



All these experimental maps show that the pressure on the (x, y) plane can be 
easily modelized by a modal decomposition of the field. On the second hand, 
for all cases, it appears that the attenuation increases along the street and 
depends on frequency : the open roof of a street canyon provides acoustic losses 
for the acoustic fled in the street. In regards to the experiment conclusions, 
two modelizations of the acoustic propagation along a street are proposed in 
the following section : an analytical model and a numerical simulation with a 
description of the attenuation due to the open roof. 
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3 Analytical modelization of the open roof effect on acoustic prop- 
agation along a street 



The analytical model of the acoustic propagation along a street is based on the 
modal decomposition of the pressure field in the horizontal plane (x, y) (see fig. 
4a). The attenuation of the pressure along a street, due to acoustic radiation 
losses through the open roof is described by means of a 2D model propagation 
in free field in the vertical plane (x, z) (see fig. 4b). The combination of these 
two approaches allows to elaborate a 3D analytical model of the acoustic 
propagation along a street. 



3.1 Acoustic propagation in a 2D waveguide 



A semi-infinite 2D wave guide of width d closed by a rigid wall at x = con- 
taininK^macoustic source is considered (see fig. 4a). In the approximation of 



R r 




R 



(a) (b) 
Fig. 4. View of the street in the (x, y) plane (a) and in the (x, z) plane (b). 



the linear acoustic and in the case of a perfect fluid, the acoustic pressure p(t) 
and the acoustic velocity v(t) satisfy the equations of the mass conservation 
law 



dp(t) 
dt 



clp V.v{t) = 0, 



and of the impulsion conservation (Euler equation) 
dv{t) 



Po- 



dt 



Vp(t) = 0, 



(2) 



where c is the sound celerity and p is the air density. In the frequency domain 
(with a temporal convention e^', u> being the acoustic pulsation), eqs. (1) and 
(2) allow to write the propagation problem, outside the acoustic source and 
with rigid boundaries, under the following form 
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(A + k 2 )p = 0, 



(3) 
(4) 




where p is here complex variable, k = uj/cq is the wave number and d/dn is 
the normal wall derivative. The solution of the acoustic problem given by eqs. 
(3) and (4) is written as 



where k 2 x = k 2 — (nir/ d) 2 , 5 n0 is the Kronecker symbol (S n0 = 6(n) = 1 for 
n = and 5 n Q = OVra ^ 0) and N is the number of modes. To determine 
the amplitude A n of each mode, the boundary condition given by the acoustic 
source is used. We suppose that the pressure field is calculated in the horizontal 
plane containing the source. The source is represented by a hole with a width d s 
in a rigid wall on the street entry and located at the point (0, y s ) belonging to 
the (x, y) plane (see fig. 4a). The source condition is described in the following 
form 



The projection of the source conditions given by eqs. (6) and (7) on the modal 
basis allows to determine the modal amplitude of each mode A n . 

3.2 Determination of the attenuation due to the open roof 

To modelize the attenuation of the pressure field due to acoustic radiation 
losses through the street open roof, a 2D propagation model is used in the 
vertical plane (x,z) (see fig. 4b). In this plane, the ground of the street is 
considered as perfectly rigid with a reflexion coefficient R g = 1 and the source 
is described by a point source located at the height h s embedded in a rigid 
wall with height h. The acoustic radiation condition describing the open roof 
of the street is modelized by means of a reflexion coefficient R r at the height 
h on the z-axis. 

The pressure field in the street can be decomposed as an infinite sum describing 
the multiple reflexions on the ground and on the street roof. In a 2D domain, 
the Green function G(f, r ) is written as 




n=l 



(5) 



v(y) = 1, for y s - d s /2 < y < y s + d s /2, 
v(y) = 0, for < y < y s - d s /2 and y s + d s /2 < y < d. 



(6) 
(7) 



G(r,r ) = - J -(H^k\f-f \)), 



8 



8 



where r*o defines the position of the source and Hq is the Hankel function of 
first order. According to eq. (8), the pressure field, at the altitude z in the 
street, takes the following form 



p(x, z)=~ Yl ( R r) lml H l ( k \/ x2 + ( z ~ 2mh ~ h s 



m=—oo 



+Hl (k^x 2 + (z + 2mh + h s )^j 



(9) 



where m informs on the number of image sources used in the model. 



3.3 Analytical model of the acoustic propagation in a street 



The complete modelization of the pressure field along the street is established 
into two steps : 

(1) the first step consists to calculate the pressure field along the street in 
the horizontal plane (x, y) without any attenuation by means of the eq. 
(5) and after having determined the modal amplitude of each mode by 
using the source conditions described by the eqs. (6) and (7); 

(2) in the second step, the attenuation of the pressure field for an altitude z, 
described by the eq. (9) is applied to each point of the pressure field in 
the horizontal plane (x,y). For this calculus, the height of the source h s 
and the altitude z of the horizontal plane in which the pressure field is 
determined are necessary. 

Finally, the pressure field along the street, at the altitude z, taking into account 
the acoustic radiation losses due to the open roof is given by the following 
equation : 



p{x,y,z) = -jJ2J2( R r) lml {A n ^2-5 n0 cos(^-y)e 
H] (k^Jx 2 + (z- 2mh - hsp) + Hi [k^x 2 + (z + 2mh + h s ) 2 ^j J (10) 



In this work, the reflexion coefficient R r describing acoustic radiation leakage 
due to the street open roof is determined by fitting the modelization with the 
experimental results obtained with the street scale model. 
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4 Numerical simulation of the open roof effect on acoustic propa- 
gation along a street 



The acoustic pressure field along the street with a height h is simulated using a 
two-dimensional finite difference time domain (FDTD) computation modified 
to take into account the acoustic radiation losses through the street open roof. 

To obtain this adaptated FDTD model, a 3D-description of the acoustic prop- 
agation is reduced to an equivalent 2D description. In this 2D model, the at- 
tenuation (called also leakage) due to the third dimension (here the z-axis) 
is introduced by means of negative source term q(t) in the mass conservation 
law (eq. (1)) leading to the following form 

^p- + clp V.v(t)=p clq(t). (11) 



This negative source, uniformly distributed in the 2D plan, depends on the 
pressure p(t) and can be written as 

q(t) = -ap(t), (12) 



where a > is a coefficient to be determined. This source term can be con- 
sidered as the effect of the third dimension (attenuation or leakage) on the 
propagation of the wave in the 2D space (x, y). 

The discretized forms of the eqs. (11), (12) and the Euler equation given by 
the eq. (2) can be achieved by a two dimensional finite difference time domain 
(FDTD) computation. After the integration of the mass conservation along a 
surface element of dimension dx along the x-axis and dy along the y-axis, we 
obtain 

^ 0, y) + CqPo f ^ 0, V) + ^ (x, y) J = -p c 2 ap(x, y), (13) 



where v x and v y are respectively the projectionss of the acoustic velocity along 
the x and y axis. The same approach is used with the Euler equation (eq. (2)) 
and a temporal discretization of these two equations allows to obtain a 2D- 
FDTD computation [19] of the propagation in a 3D space. 

To apply the adaptated 2D-FDTD simulation to the propagation along a 
street with an open roof, the discretized equations are computed by means 
of a Matlab program where the boundary conditions are introduced on the 
acoustic velocity (the facades of the street are considered as perfectly rigid) 
and the source condition is introduced on the pressure. The coefficient a in 
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the equation (13) can be considered as an attenuation for the 2D propagation 
in the (x, y) plane and describes the acoustic leakage due to the open roof of 
the street. This coefficient depends on the altitude z of the 2D plan (x,y) in 
comparison with the source height h s in the street, on the aspect ratio of the 
street d/h and on the acoustic radiation condition on the street roof. 



5 Results and discussion 

5.1 General observations 

In this section, we propose a qualitative comparison between the experimental, 
analytical and numerical results for / = 1000 Hz, / = 1500 Hz and / = 2500 
Hz. 

Two simulated results of the acoustic pressure field in a street with the 
adapted-2D FDTD computation are proposed. The simulation modelizes the 
experimental source condition and a Perfect Matching Layers (PML) is used 
to describe the anechoic terminaison of the street. The comparison with the 
experimental results is led here for / = 1000 Hz and / = 2500 Hz. The coef- 
ficient a modelizing the acoustic losses due to the street open roof is adapted 
by means of a minimization algorithm with the experimental data. 

The results of the analytical model are proposed for / = 1000 Hz and / = 
1500 Hz. The acoustic radiation conditions of the open roof are described 
by a reflexion coefficient in the model. As planned, the reflexion coefficient 
is negative and its value is adapted by fitting the analytical results with the 
experimental ones. The analytical field is calculated by means of a modal 
decomposition with N modes and the attenuation along the street is modelized 
with m + 2 image sources (see eq. (9)). 

The fig. 5b proposes the analytical pressure map in the street computed by 
means of the eq. (10) with N = 56, m = 3 and R r = —0.3 for / = 1000 Hz 
with a square source centered on y s = 0.175 m and at the altitude z = 0.07 
m. The qualitative agreement with the corresponding experimental result (see 
fig. 5a) is correct : the analytical shape of the field is close to the experimental 
one. For this case, at the end of the street, the location on the y-axis of the 
minimum of the analytical map are shifted refering to the experimental results. 
This difference can be explained by a default in the alignement of the robot 
and the street x-axis. 

Fig. 5c shows the simulation of the acoustic pressure along the street for 
/ = 1000 Hz with a square source centered on y s = 0.175 m. The step of 
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spatial sampling of the FDTD is 0.01 m. The comparison between fig. 5c 
and 5a shows a good agreement between both results. The simulated result 
exhibits the same shape of the acoustic field than the experimental one. The 
attenuation of the simulated and experimental results are in the same order 
which prooves that, for this frequency, the acoustic radiation conditions for 
the street roof can be modelized by an attenuation coefficient. 
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Fig. 5. (a) Experimental acoustic pressure in a street (scale model) for / = 1000 
Hz at z = 0.07 m and for a square source centered on y s = 0.175 m and z s = 0.07 
m, (b) analytical acoustic pressure in a street at z = 0.07 m, for / = 1000 Hz and 
for a square source centered on y s = 0.175 m., (c) adaptated 2D-FDTD simulation 
of the acoustic pressure in the 2D-street for / = 1000 Hz and for a square source 
centered on y s = 0.175 m. 

The fig. 6b shows the analytical calculation of the acoustic pressure in the 
street at z — 0.07 m for a frequency / = 1500 Hz with the same source 
condition than in the experimental study. The number of modes is N = 59, 
the value of the reflexion coefficient is estimated to R r = —0.1 and we use 5 
image sources. The qualitative comparison of this result with the experimental 
one (see fig. 6a) shows a good agreement except at the end where the analytical 
model overestimates pressure compared to experiments. 




1 / Y 1.5 

x (m) 

Fig. 6. (a) Acoustic pressure in the street (scale model) for / = 1500 Hz at z = 0.07 
m and for a square source centered on y s = 0.175 m and z s = 0.07 m, (b) analytical 
pressure map in a street at z = 0.07 m for / = 1500 Hz and for a square source 
centered on y s = 0.175 m. 



The fig. 7 proposes the simulated map of the acoustic pressure along the street 
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for / = 2500 Hz. The step of spatial sampling of the FDTD is 0.01 m. The 
comparison of this result with the experimental map shows a good agreement. 
The mode distribution of the pressure field along the street is the same for 
the experimental and simulated results. The acoustic radiation conditions are 
taken into account with qualitatively good precision. 




Fig. 7. (a) Acoustic pressure in the street (scale model) for / = 2500 Hz at z = 0.07 
cm and for a square source centered on y s = 0.1 m and z s = 0.07, (b) adaptated 
2D-FDTD simulation of the acoustic pressure in the 2D-street for / = 2500 Hz and 
for a point source centered on y s = 0.1 m. 

In view of these primary results, the effect of the street open roof can be 
modelized by a attenuation coefficient in the 2D adapted FDTD computation 
depending on the frequency : for / = 1000 Hz, a = 0.0018 s -1 and for / = 2500 
Hz, a = 0.0025 s" 1 . The attenuation is greater at 2500 Hz than at 1000 Hz. For 
the analytical model, the open roof effect is described by a negative reflexion 
coefficient with a amplitude which decreases when the frequency increases 
: Rr = -0.3 for / = 1000 Hz, R r = -0.1 for / = 1500 Hz and R r = for 
/ = 2000 Hz and 2500 Hz. This result shows that more the frequency increases, 
more the open roof radiation conditions are close to a complete absorption. 



5.2 Quantitative comparison : low frequency validation 



To study with more precision the differences between the experimental, mod- 
elized and simulated results, the acoustic pressure along the street for y = 0.09 
m and y = 0.15 m is compared for / = 1000 Hz on figs. 8a and 8b. The analyt- 
ical model proposes a good prediction of the pressure attenuation and of the 
maximum pressure positions along the street. With the adaptated 2D-FDTD 
simulation, the attenuation is well predicted and the maxima positions of the 
pressure field along the street are estimated with a maximum error of 0.05 m. 
This result prooves that a reflexion coefficient used in the analytical model 
and a negative pressure source uniformly distributed along the street used 
in the 2D-FDTD simulation constitute a good modelization for the pressure 
attenuation along the street due to the open roof. 

The experimental and analytical acoustic pressure profiles along the street for 
/ = 1500 Hz at y = 0.05 m are proposed in fig. 9. This comparison shows 
a good agreement between the analytical model and the experimental study. 
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Fig. 8. Comparison of the experimental (blue curve), analytical (black curve) and 
numerical (gray curve) results of the pressure profile along the street for / = 1000 
Hz at y = 0.09 m (a) and for / = 1000 Hz at y = 0.15 m (b). 

The maximum locations of the pressure along the street are well estimated 
and the global attenuation is predicted with a good accuracy. This prooves 
that the use of a negative reflexion coefficient can be a good modelization of 
the acoustic radiation conditions by the open roof at this frequency. 
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Fig. 9. Comparison of the experimental (blue curve) and analytical (black curve) 
results of the pressure profile along the street for / = 1500 Hz at y = 0.05 m. 



5.3 Frequency limit of the modelization 



For / = 2000 Hz, the analytical, simulated and experimental results are com- 
pared on the fig. 10 for y = 0.08 m. The analytical model is used with R r = 
and N = 62 and the simulated one is performed with a = 0.0020. The agree- 
ment between the simulated and experimental curves is good : the attenuation 
is very well predicted and a error of 0.03 m on the maximum positions along 
the street, involving notably by the discretization of the FDTD computation, 
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is visible. 



For the analytical results at / = 2000 Hz, the limits of the model are reached. 
The positions of the maximum are well estimated with this method but the 
attenuation is badly predicted. The reflexion coefficient is chosen to zero which 
involves a maximum for the pressure attenuation. But the comparison with 
the experimental data shows that the analytical attenuation is underestimated 
: at x = 2.5 m the error between the model and the experience for the acoustic 
pressure is greater than 3 dB. 




X (m) 

Fig. 10. Comparison of the experimental (blue curve), analytical (black curve) and 
numerical (gray curve) results for the pressure attenuation along the street for 
/ = 2000 Hz at y = 0.08 m. 

The fig. 11 shows the experimental, analytical and simulated acoustic pressure 
along the street for / = 2500 Hz at y — 0.1 m. The simulated result is obtained 
with a = 0.025 and the analytical model is performed for R r = with 11 image 
sources (corresponding to m = 5). As for / = 2000 Hz, the agreement between 
the numerical and experimental results is good : the pressure attenuation 
and the maximum positions of the pressure field are well predicted excepted 
the near foeld of the source. The analytical model overestimates the pressure 
attenuation along the street but proposes a good estimation of the maximum 
positions. 

The analytical model or the numerical simulation, developped in this work, 
are based on a low frequency hypothesys which permits to modelize the at- 
tenuation by a single coefficient depending on the frequency and the street 
size. By increasing the frequency, the acoustic radiation losses, due to the 
open roof of the street, is not anymore uniformly distributed on all modes 
in the waveguide. The notion of leaky modes [20] (with its own attenuation 
or leakage) must be introduced in the modelization to describe the acoustic 
propagation along the street. The fig. 12 shows the map of the experimental 
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Fig. 11. Comparison of the experimental (blue curve), analytical (black curve) and 
numerical (gray curve) results of the pressure profile along the street for / = 2500 
Hz at y = 0.1 m. 

pressure amplitude along the street for a frequency / = 3400 Hz. 
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Fig. 12. Experimental acoustic pressure in a street at z = 0.07 m, for a frequency 
/ = 3400 Hz and for a square source centered on y s = 0.1 m and z s = 0.07 m. 

On the fig. 13, the pressure profile for / = 3400 Hz at y = 0.1 m is represented. 
The difference of leakage between the modes is easily shown at the beginning 
of the street : after 1 m along the street, a leaky mode (defined by 2 A = d 
along the y-axis) disappears completely with a global attenuation illustrated 
by the shaded line while a second leaky mode (defined by A/2 = d along the 
y-axis) propagates with a different attenuation illustrated by the second line. 



6 Conclusion 

In this work, we have shown that it is possible to modelize and to simulate the 
acoustic propagation along an urban street with a good agreement comparing 
to experimental results. For the low frequency case (until 100 Hz at full scale), 
an adaptated 2D-FDTD simulation, taking into account the acoustic losses 
due to the radiations by the street roof, proposes a good modelization of the 
acoustic propagation. A coefficient a of the attenuation along the street is 
adapted by fitting with the experimental results. 

An analytical model for the acoustic propagation along a street canyon is also 
developped by mixing a modal approach to describe the propagation in the 
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Fig. 13. Experimental pressure profile along the street for / = 3400 Hz at y = 0.1 
m. 

horizontal plane and a 2D propagation model using image sources to modelize 
the attenuation along the street. A reflexion coefficient is defined to modelize 
the role of the street open roof. For the low frequency case, the decrease of the 
reflexion coefficient modulus with the increase of the frequency shows that the 
acoustic radiation conditions of the open roof depend on frequency and prooves 
that more frequency increases less reflexions on the open roof are present. For 
higher frequencies, while the adaptated 2D-FDTD simulation describes well 
the pressure field, the analytical model is limited by the presence of numerous 
leaky modes with different attenuations. 
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